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ABSTRACT 



If the observed light neutrino masses are induced by their Yukawa couplings to singlet right- 
handed neutrinos, natural smallness of those renders direct collider tests of the electroweak scale 
neutrino mass mechanisms almost impossible both in the case of Dirac and Majorana (seesaw of 
type I) neutrinos. However, in the triplet Higgs seesaw scenario the smallness of light neutrino 
masses may come from the smallness of B — L breaking parameters, allowing sizable Yukawa 
couplings even for a TeV scale triplet. We show that, in this scenario, measuring the branching 
fractions of doubly charged Higgs to different same-charged lepton flavours at LHC and/or ILC 
experiments will allow one to measure the neutrino mass parameters which neutrino oscillation 
experiments are insensitive to, including the neutrino mass hierarchy, lightest neutrino mass 
and Majorana phases. 
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1 Introduction 



In recent past neutrino oscillation experiments have shown convincingly that at least two light 
neutrinos have non-zero masses and their mixing is characterized by two large mixing angles 
PQ. Those facts constitute indisputable evidence of new physics beyond the standard model 
(SM). However, despite of the intense experimental and theoretical effort over many years, 
understanding of the origin of neutrino masses is still missing. 

From the experimental point of view the information on neutrino masses coming from oscil- 
lation experiments is limited by the fact that these experiments are only able to measure the 
differences of squared neutrino masses and not their absolute magnitude, neither are they sen- 
sitive to the Dirac or Majorana nature of light neutrinos. In particular, the present oscillation 
experiments cannot distinguish between the two possible mass ordering patterns of light neu- 
trinos, the normal and the inverted ones, and are insensitive to the possible Majorana phases 
[2] of neutrinos. The observed smallness of neutrino #13 mixing angle makes it very difficult 
to measure any new parameter, such as the neutrino Dirac CP phase 5, in neutrino oscillation 
experiments before a distant-future neutrino factory [3]. To learn conceptually new facts about 
light neutrinos in a shorter time-scale requires likely an experimental breakthrough either in 
low energy neutrino experiments, such as the neutrinoless double beta decay (OuP/3) decay 
experiments, or in collider physics. 

From the theory side we still do not know why neutrinos are so light compared to charged 
fermions. It is natural that the 577(2) l doublet neutrinos couple to new singlet (right-handed) 
neutrinos N and the SM Higgs doublet in a direct analogy with all other SM fermions. If 
this is the only new physics, the smallness of neutrino masses requires unnaturally small Dirac 
Yukawa couplings. Alternatively, the singlet neutrinos may have very large Majorana masses 
which suppress the light neutrino masses to the observed range via the seesaw mechanism of 
type I [U El El E] even for large values of the Yukawa couplings. Generically neither of those 
simple scenarios can be directly probed at low energy nor collider experiments^. Making the 
singlet neutrinos as light as 1 TeV to be kinematically accessible at colliders does not help 
because their only interactions are of Yukawa type and the seesaw mechanism predicts that 
the couplings are too small for any observable signal except the neutrino masses. Complicated 
model building is required to ensure the correct light neutrino masses, 1 TeV heavy neutrinos 
and meaningfully large neutrino Yukawa couplings at the same time. Unfortunately the direct 
tests of singlet neutrino mass mechanism at LHC are experimentally demanding even in those 
models [El [131 E] . 

However, group theory tells us that generation of non-zero masses for the SM doublet neutrinos 
does not require the existence of singlets. One of the best motivated and best studied neutrino 
mass scenario is the triplet Higgs mechanism [15], [El El EH1 EE], sometimes called seesaw 
mechanism of type II. From the point of view of direct tests the triplet neutrino mass mechanism 

1 This conclusion may be different if softly broken supersymmetry exists in Nature. Flavour violating Yukawa 
couplings of heavy neutrinos may induce flavour off-diagonal elements in the soft slepton mass matrices via the 
renormalization effects [TU] which may lead to observable rates of lepton flavour violating processes. This 
very complex scenario requires analyses beyond the present one [llj . 
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has several advantages over the singlet one. Firstly, the SU(2) L triplet multiplet contains a 
doubly charged scalar which can be pair produced at colliders independently of their Yukawa 
couplings. Thus tests of this mechanism are limited only by the collision energy. Secondly, 
the smallness of neutrino masses does not imply the smallness of triplet Yukawa couplings. 
As neutrino masses in this scenario are necessarily of Majorana type, they may be different 
from the Dirac fermion masses because of the smallness of B — L breaking. This is natural 
by the 't Hooft criterion as B — L is a conserved quantum number in the SM. Thus the 
neutrino Yukawa couplings to triplet may be sizable, constrained by unobserved lepton flavour 
violating interactions, and dominate over the triplet coupling to two gauge bosons. Thirdly, 
the triplet Yukawa couplings directly induce the neutrino mass matrix up to the small B — 
L breaking triplet vacuum expectation value (VEV) which appears in neutrino masses as a 
common proportionality factor. Altogether those arguments imply that one can study the 
neutrino mass parameters at Large Hadron Collides (LHC) and/or International Linear Collider 
(ILC) experiments by just counting flavours of the same-charged lepton pairs originating from 
the doubly charged Higgs boson decays. 

In this work we extend the analyses of our previous paper [2D]- While in Ref. [20J we studied 
the discovery potential of LHC experiments for the process pp — > [21] [22], [23] EU [25], 

[26] [20] [27] assuming that the subsequent decays of $ ±± are determined by neutrino data, in this 
analyses we turn the argument around and study what can one learn about neutrino physics 
if LHC and/or ILC will discover the triplet Higgs bosons. In particular, we concentrate on 
neutrino parameters which cannot be measured in oscillation experiments, the light neutrino 
mass ordering, the mass of the lightest neutrino and the Majorana phases ai and a<i- Neutrino 
mass hierarchy patterns at colliders have been previously studied in [28]. First we derive 
analytical expressions for those quantities which are functions of the doubly charged Higgs 
branching fractions to different flavour combinations of charged lepton pairs, $ ±± — > £i£j, 
i,j = e,fi,r. Thus neutrino physics at colliders turns out to be just a counting experiment of 
lepton flavours. This simplifies the life in particular at LHC experiments which, in general, have 
larger measurement errors than at ILC. The analytical results are first derived assuming the 
tri-bi-maximal mixing for neutrinos, which predicts sin #13 = 0, and extended later to non-zero 
values of sin# 13 . After that we study to what precision those quantities can be measured in 
realistic experiments. Finally we demonstrate that combining positive colliders signals of this 
scenario with the possible measurement of the neutrino mass matrix entry (m u ) ee would allow 
one to determine separately the size of triplet Yukawa couplings and the B — L breaking VEV of 
the triplet. Thus one can entirely probe the neutrino mass generating mechanism at terrestrial 
experiments. 

We find that there are distinctive flavour signals which indicate certain patterns of neutrino 
mass matrix. For example, very few electrons in $ decays definitely points towards normally 
hierarchical light neutrinos. There is theoretical ambiguity in determination of the Majorana 
phases and only a combinations of them can be measured. However, in the case of very hi- 
erarchical neutrino mass spectrum one of the Majorana phases is effectively decoupled from 
physics and one can, in principal, measure the magnitude of the physical phase. Although the 
experimental errors in determining those quantities may turn out to be quite large in general, 
we show that there exist scenarios which can already be fully solved at LHC. In the optimistic 
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scenarios the branching fractions of the doubly charged Higgs boson decays can be used to (z) 
determine the neutrino mass hierarchy; (ii) estimate the mass of the lowest neutrino state; (Hi) 
estimate the Majorana phases of CP violation; (iiii) measure the value of Higgs triplet VEV. 
We note that those measurements are also sensitive to all other neutrino parameters including 
the mixing angles and CP violating phase 5. We show that the latter two are, in principle, 
measurable at collider experiments. However, those quantities can be determined with much 
higher precision in other experiments and we do not study their effects in detail in this paper. 

The paper is organized as follows. In Section[2]we present the collider phenomenology of doubly 
charged Higgs boson and relate the collider observables to the neutrino parameters. In Section 
[3] we present details of the analysis of neutrino parameter measurements at colliders. In Section 
H] we discuss the possibility of measuring the triplet Higgs VEV and determination of the full 
neutrino mass matrix. Finally we conclude in Section [51 



2 Phenomenological setup 

In this work we assume that the SM particle spectrum is extended by a scalar multiplet $ 
with the SU(2)l x U(1)y quantum numbers $ ~ (3, 2). We also assume that its mass is below 
0(1) TeV and the pair production processes at colliders, 

pp —> $++$ and e + e _ — > $ ++ $ , (1) 

are kinematically allowed. Such a scenario is realized, for example, in the little Higgs models 

I2H EOl EH E2] - 

The triplet couples to leptons via the Lagrangian 

L = i£ D Li T 2 Yl j (T-$)l Lj + h.c., (2) 

where are the Majorana Yukawa couplings of the triplet to the lepton generations i,j = 

e,fi,r. If the neutral component of triplet acquires a VEV v$, the non-zero neutrino mass 
matrix is generated via 

K)jj = 2(Ya f ) ij Vis>. (3) 

To avoid the existence of phenomenologically unacceptable Majoron the B — L breaking VEV 
v§ cannot occur spontaneously. Instead it should be induced effectively via the coupling of $ 
to the SM Higgs doublet H as fi^°H°H°, where the dimensionful parameter /i breaks B — L 
explicitly [33] • Because in the limit \i — > the symmetry of the model is enhanced, it is 
natural that fi is a small parameter. Model building in this direction [Ml 1351 136] is beyond 
the scope of the present analyses. Indeed, the above described scenario is consistent with the 
observation that neutrino masses are much smaller than the masses of other SM fermions. Thus 
the smallness of neutrino masses is explained by be smallness of v$ and the Yukawa couplings 
(Y$)ij can be of order SM Yukawa couplings. The most stringent constraint on them arises 
from non-observation of the muon decay \i — > eee which implies Y ee Y*^ < 2 • 10~ 5 for m$ = 1 
TeV [37]. 
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It is important to emphasize that the precise values of are not relevant for the collider 

physics we consider in this work. The relationship between neutrino parameters and doubly 
charged Higgs boson decays comes from the fact that the Yukawa coupling matrix of doubly 
charged Higgs to leptons is proportional to the Majorana mass matrix as given by Eq.([3]). Thus, 
to establish this connection experimentally, observable rates of the leptonic branching fractions 
must exist. 

The decay width of doubly charged Higgs to the corresponding leptonic channel is given by 



k\(Y$)ii\ 2 mg,±± i = j, 



r y = r($ ±± - ifif) = { (4) 

and the decay width to the WW channel is 

iw - r ( $- wtwt) = ags^ ( 3 <t + ^- l ) ( 1 .1<£f s kvl 

167im w ± \ m U± 4m w ± J V m * ±± / 

(5) 

The branching ratio of $ ±=l= to a single leptonic channel can be calculated using the decay 
widths 

BR,, = BR($ ±± - £f£f) = (6) 

1 tot 

where Y tot = J2i>j + ^ww is the total decay width. Since T^- is directly related to neutrino 
mass matrix, we can derive a relation between the $ ++ $ branching ratios that can be 
measured in collider experiments and the neutrino mass matrix, that contains all currently 
unknown neutrino parameters. The branching ratio to a single decay channel can be found 
combining Eqs. O), © and (jUJ) 



BR($ ±± - iflf) = _ J (m ^J - 4 , (7) 



where (m v )ij is the neutrino mass matrix in flavor basis. Eq. ((7|) shows the direct relationship 
between neutrino parameters and the $ branching ratios. 

The neutrino mass matrix can be diagonalized by unitary leptonic mixing matrix U, 

m v = U*m°U\ (8) 
where the diagonalized neutrino mass matrix is given by 



Here mi, m,2 and represent the masses of neutrino mass eigenstates v%, and v^. v\ and V2 
masses differ by Am^, = 7.92(1 ± 0.09) x 10~ 5 eV 2 measured by solar oscillation experiments 
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[38J and m x < m 2 . The third eigenstate z/ 3 is separated from the first two by splitting Am^ (m = 
2.6(llo'i5) x 10 _3 eV 2 [3U] and can be heavier or lighter than the solar pair. The two possibilities 
are called normal and inverted spectrum, respectively. The third possibility - nearly degenerate 
masses - appears when the lowest neutrino mass is large compared to the measured mass 
differences and mi « m 2 ~ m 3 . Cosmology implies that neutrinos are lighter than about 0.2 
eV [10]. 

Since we have assumed that there are only three Majorana neutrinos, U is a 3 x 3 mixing matrix 
that depends on three mixing angles and three phases and can be parameterized as 

/I \ / c 13 s 13 e" t5 \ ( c 12 s 12 \ / e*°* \ 
U = c 23 S23 1 -s 12 c 12 e ia * , (9) 

V -S23 c 23 / V ~ s rse i5 c 13 / \ l/ \ 1 J 

where = cos 6^, = sin^- and denote the mixing angles. The quantities 5, a± and 
ct 2 are CP violating phases. 5 is the Dirac phase and characterizes CP violation regardless 
of the character of neutrinos. a± and a 2 are called Majorana phases and are physical only if 
neutrinos are Majorana particles. If the neutrinos were Dirac fermions, both Majorana phases 
could be absorbed by appropriately redefining the neutrino fields, and the only observable CP 
violation parameter would be the Dirac phase 5. Also note that 5 appears in the mixing matrix 
only as sin 9i 3 e lS - so the influence of 5 crucially depends on the value of 6*i 3 and has physical 
consequences only if 813 is non-zero. 

The mixing matrix contains six independent parameters: three mixing angles (#13, # 23 and #i 2 ) 
and three phases (a\, a 2 and 5). Mixing angles are known from the global fit to neutrino 
oscillation data and are given by (2a errors) [381 EH] 

sin 2 6 12 = 0.314(1^) , sin 2 fl 23 = 0.45(1^) , sin 2 9 13 = 0.8^-3 x 10 -a_ (10) 
Up to now no experiment has been able to determine the values of phases so that 

5, an, a 2 e [0,2vr], (11) 

remain unconstrained. 

So far we have shown that doubly charged Higgs boson decay statistics is directly related 
to neutrino parameters as given by Eq. (13) and depends on neutrino mass matrix and 
mixing matrix U. Additional information can be acquired from the neutrinoless double beta 
decay experiments, that independently probe the absolute value of the (m u ) ee entry in neutrino 
Majorana mass matrix. Such relations allow direct measurements of neutrino parameters in 
particle collider experiments. 

3 Measuring neutrino parameters at colliders 

Doubly charged Higgs boson has 6 different leptonic decay channels. Branching ratios to these 
channels are functions of neutrino parameters according to Eq. We have fixed the values of 
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mass differences Am 2 tm and Am^, in the subsequent calculations, as they have been measured 
with a good precision in neutrino oscillation experiments. With such an assumption we can 
write an equation system of six independent equations that relates branching ratios of six 
different $ ±=t leptonic decay channels with unknown neutrino parameters, 

BR„ = fk{m , sign(Am atm ), 9 13 , 6> 23 , 9 12 , 5, ax, a 2 ), (12) 

where mo represents the mass of the lowest neutrino mass eigenstate (mi or m 3 for normal 
or inverted mass spectrum, respectively), k = 1,...,6 and i,j = e,fj,,r. In the subsequent 
analyses we have used relations between the leptonic branching ratios instead of their absolute 
values. This method is independent of the possible $ ±=t decay to WW channel, which is more 
complicated to measure accurately at LHC. In such an approach we simply count the events of 
$ ±± decays to different channels and calculate their relative differences. As a result we have 
five independent equations. In order to solve this equation system with respect to unknown 
neutrino parameters, we have to fix at least some of them. Consequently we can solve the 
equation system (1121) for different neutrino parameters and obtain them as functions of the 
$ ±± leptonic branching ratios BRy. 

3.1 Results for the tri-bi- maximal mixing 

Since approximate values of neutrino mixing angles are known from oscillation experiments 
and the precision of measurements is expected to be increased in upcoming years [41J, we fix 
their values in most of our analyses. We have chosen to follow the tri-bi-maximal model [42J 
. It has been proposed that the combined existing data from neutrino oscillations point to 
a specific form of the lepton mixing matrix with effective bi-maximal mixing of and v T at 
the atmospheric scale and effective tri-maximal mixing at the solar scale - hence denoted as 
tri-bi-maximal mixing. The tri-bi-maximal mixing predicts 

sin 2 6»i2 = - , sin 2 6» 23 = ~ , sin 2 13 = 0, (13) 

which are perfectly compatible with the present experimental uncertainties given by Eq. (fTUI) . 
The main aim of this paper is to provide information about the Majorana phases and absolute 
values of neutrino masses. In the tri-bi-maximal model the CP violating phase 5 is not physical 
due to the zero value of 9\% and the only remaining unknown variables are the lowest neutrino 
mass mo, neutrino hierarchy i.e. sign(Am a j m ) and Majorana phases a± and a 2 - 

Having fixed the mixing angles according to Eq. ([TBI , we end up with four independent equations 
for branching ratios, since BR eA1 = BR er and BR W = BR TT . If a measurement would show that 
these branching ratios are not equal, this is a clear indication that the tri-bi-maximal model has 
to be modified. As we are using the relations between branching ratios for the calculations, the 
number of independent equations is reduced to three. Such equation system can be solved with 
respect to three unknown parameters: the lowest neutrino mass mo and Majorana phases a.\ 
and «2- We show how the mass of the lowest neutrino mass eigenstate, neutrino mass hierarchy 
and the difference of two Majorana phases |Aa| can be uniquely determined from the relation 
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(j7|). Unique solutions for ot\ and a 2 are not determined by the $ ++ $ branching ratios, and 
two sets of degenerate solutions are found. 



3.1.1 Neutrino hierarchy and the lowest neutrino mass 

First we consider the equation system given by Eq. (fT2j) with the fixed tri-bi-maximal mixing 
angles. For the neutrino mass hierarchy and lowest neutrino mass determination we combine 
the branching ratios of /i/z, fir, ee and e/x channels. After some simple algebra we find a relation 
between these branching ratios that depends only on neutrino masses and is independent of 
the Majorana phases, 

= 2BR W + Bfy - BR ee _ -m? + m\ + 3m 2 3 
1 ~ BR ee + BR eM 2ml + m 2 2 ' [ ' 

Here and onwards in this paper C x denote constant dimensionless parameters which can be 
measured in experiments. 

The mass hierarchy can be easily determined by simply measuring the value of C\ that is 
independent of the values of aci and a.%. It can be found that G\ uniquely determines the mass 
hierarchy as follows: 

• C\ > 1 - normal mass hierarchy, 

• Ci < 1 - inverted mass hierarchy, 

• G\ ~ 1 - degenerate masses. 

After the mass hierarchy measurement we can solve Eq. (TH|) for either normal or inverted mass 
hierarchy. For the normal mass hierarchy m\ is the lowest mass state, m\ = ml + Am^ and 
m\ = ml + Am^ ( + Am^ lm . After substituting m 2 and m 3 , we get the following equation that 
can be solved with respect to mi, 



AmL(4-C 1 ) + 3Am, 



at in 



3(C! - I) 



(15) 



a t in 



Alternatively, for the inverted mass hierarchy m 3 is the lowest mass state, mi, = m 2 , + Am 
and ml = m\ + Am^ m — Am 2 ol . After the substitutions, Eq. ( THl) can be solved with respect 
to m^ as follows: 

2 Amg oi (l + 2Ci)-3CiAmg tm 

3(d-l) ' [ ' 

For nearly degenerate masses (mi > O.leV) accurate measurement of the lowest neutrino 
mass requires very good experimental precision (which is not likely to be achieved at LHC) 
because the branching ratios are increasingly less mass dependent for larger mass values. This 
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Branching ratios (normal hierarchy) Branching ratios (inverse hierarchy) 




Figure 1: Distribution of the $ ++ leptonic branching ratios as a function of the lightest neutrino 
mass. The left (right) panel corresponds to the normal (inverted) mass hierarchy. For nearly 
degenerate masses the two possibilities imply almost the same result. 



is demonstrated in Figure [T] which presents the dependency of doubly charged Higgs branching 
ratios on the lightest neutrino mass for the normal and inverted mass hierarchies. We have 
assumed a real mixing matrix i.e. fixed Majorana phases to zero. The e/x and er channels have 
only vanishingly small contributions for a± = ct 2 = 0, but are increased for non-zero values of 
the Majorana phases. The branching ratio to ee channel is a especially good characteristic for 
mass hierarchy determination that varies greatly depending on the hierarchy and the neutrino 
mass. This branching ratio is negligible for the normal mass hierarchy with very small mass 
while it is the dominant decay channel for the inverted mass hierarchy. If the mass of the lightest 
state increases, both the normal and inverted hierarchies have almost the same distribution of 
branching ratios, $ ±=t decay to ee, /i/i and rr with nearly equal probabilities while the decays 
to other channels are negligible. This indicates the degenerate masses. 



3.1.2 Majorana phases 

If the neutrino masses are measured as shown in the previous section, we can determine the 
values of Majorana phases in a similar way. Once again we use the tri-bi-maximal values for 
all mixing angles and combine expressions from the equation system ( Fl2l) . We first determine 
the difference between the Majorana phases Act = [ati — CK2I- Using a relation between the ee 
and e/i decays channels we obtain 

^ _ BR efl _ 2(m\ + m 2 2 - 2mim 2 cos Act) 
BR ee Am\ + m\ + 4mim2 cos Act 

From this expression we can find separate solutions for the different mass hierarchies. For the 
normal hierarchy Act can be found to be 



(4 - 5C 2 )m2 + (2 - C 2 )Am 



2 

sol 



Act = arccos ^ zy 1 v , zy s oi , (1 

1 4(l + C 2 )m lv /m? + Am^ 1 
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while for the inverted hierarchy we find 

Aa = arccos ( 2(2C 2 1)A^ L + (4 - 5C 2 )(A< m + ^|) 

1 4(1 + C 2 )y/(Ami tm + mi)(mi + Am^ - Am'jJ 1 



Aa = arccos — -f- + 0[ f . (20) 



which can be approximated as 

4~5C 2 \ Q /Ang 

For the inverted hierarchy, up to small corrections, the equation for Aa is independent of the 
value of m 3 . This means that Eq. (|20|) is valid both for the inverted hierarchy and degenerate 
mass spectrum. The solution for the normal hierarchy given by Eq. fll8p contains the lowest 
neutrino mass which must be measured previously with an acceptable precision. 

We found that |Aa| can be uniquely determined up to a sign uncertainty sgn(ai — a 2 ) since 
cosine is an even function. In order to find a solution that separately determines ot\ and a 2 , 
we use the expression for Aa given either by Eq. (|TH|) or Eq. (I2"01) together with the definition 
of cosine of the difference of angles and construct the equation system of two independent 
equations, 

_ 261-1^ — 611^,- 2777,3(00801/71! + 2 cos a 2 m 2 ) 

' 3 



BR ee + BR e/i 2m{ + m 2 2 

cosAa = cosai cosa 2 + sinai sina 2 . (21) 

Unfortunately such an equation system does not have a unique solution due to the uncertainty 
in sgn(ai — a 2 ) and two sets of degenerate solutions for ai and a 2 are found, one of which is 
correct for a\ > a 2 and the other corresponds to a 2 > ol\. It is not possible to tell only from 
the collider data which angle is bigger and which of the solutions is correct. In the following we 
present the effect of the Majorana phases to branching ratios for three different mass hierarchies 
and discuss the consequences of non- vanishing # 13 . 



3.2 Measuring Majorana phases for different mass hierarchies 

In this section we study some well motivated particular cases of neutrino mass parameters 
which can be well measured at LHC Those are: 

• Normal mass hierarchy, m,\ = 0. Branching ratios are independent of a^ a 2 can be 
determined. 

• Inverted mass hierarchy, 7773 = 0. Branching ratios are independent of absolute values of 
Majorana phases; Aa = \a% — a 2 | can be determined. 

• Nearly degenerate masses, nil w m 2 pa m 3 = m. Expressions for branching ratios become 
independent of m. 

We have kept the mixing angles fixed to tri-bi-maximal values unless stated otherwise. 
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Figure 2: Distributions of the branching ratios as a function of a 2 - The left panel corresponds 
to #13 = with the ee, e/i and er channels giving nearly negligible contributions. When #13 
is non-zero (the right panel), small branching ratios to e/i and er channels can be measured. 
Non-zero S in the right panel causes the slight asymmetry with respect to a 2 = 7T- 



3.2.1 Normal hierarchy, mi = 



In this case the doubly charged Higgs branching ratios are independent of a\ and we can 
determine a 2 . Branching ratios to the decay channels that involve electrons can be neglected 
and, for expressing the solutions, we use the branching ratios to \x\i and fir channels. Relation 
between these channels gives the following equation with a 2 as the only unknown parameter, 



C A = BR ^ = 13Am L + 9Am lt m + 12 cosa 2 ^Am 2 sol ^Am 2 sol + AmJ^ 
BR MT 2(13Am^ + 9Am 2 atm - 12 cos a 2 ^ Am 2 sol ^ Am 2 sol + Am^J ' 

This can be solved uniquely for a 2 as 

-. I (2^-l)(13AmL + 9A<J 

12(l+2C 1 ) x /A^ v /Am^ + AmJ tm 



The distribution of branching ratios for the tri-bi-maximal mixing angles is shown in the left 
panel of Figure [21 The dominant decay channels are /z/i, rr and /it. Decays including electrons 
can be neglected, since the branching ratios to the ee, e/i and er channels are suppressed by 
Am L«/ Am atm t na t is small compared to the relevant terms in other decay channels. Non-zero 
ot2 causes a small variation in branching ratios, the fir channel is increased for while the rr 
and fifi channels are reduced proportionally. The right panel shows the effect of non-zero #13 
and 5 that create small non-zero contributions to the e/i and er channels. In this case non-zero 
6^3 could be clearly detected. However, those can comprise only about 10% of all the decays 
and require high statistics to be adequately measured at colliders. We also emphasize that the 
asymmetry of distributions in this case is CP-violation effect due to non- vanishing Dirac phase 
S. 
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In conclusion, if we have identified the normal mass hierarchy with nearly zero value of mj, 
which is being described by /i/i, fir and rr as the dominant decay channels, we can measure a<i 
from the ratio between \x\x and \it channels. We note that the changes in branching ratios are 
symmetrical with respect to 0:2 = ^ and we always have two possible solutions. The non-zero 
#13 and 5 can create a slight asymmetry in the solutions due to the CP-violation and thus 
provide a possibility of unique determination of a^- However, this is a very small effect that 
requires a precision measurement and most likely can not be detected at the LHC. 

3.2.2 Inverted hierarchy, TO3 = 

In this case the branching ratios do not depend on absolute values of Majorana phases and 
only their relative difference Aa can be measured. We still can use Eq. (1201) to determine the 
value of Aa. The distribution of all $ ±± branching ratios as functions of Aa is presented in 
Figure [3j One can see that the changes in branching ratios caused by non-zero Aa are much 
more prominent than the changes caused by a<i for normal hierarchy. When Aa = 0, ee is the 
dominant decay channel and the e/x and er channels can have only very small contributions 
resulting from small non-zero #13. If #13 = is assumed, any non-zero contribution to the e/i 
or er channels would indicate non-zero value of Aa. Non-zero Aa suppresses the ee channel 
considerably and the branching ratios to e/i and er channels can occupy more than 80% of all 
leptonic decays. Branching ratio to the ee channel remains non-zero in this case. The effect of 
non-zero #13 and CP violation angle 5 is presented in the right panel of Figure El It makes the 
distributions slightly asymmetric with respect to Aa = 7r which, in principle, can be measured. 
As in the previous case, the asymmetry is a signal of CP-violation. 



Branching ratios (8 13 = 0) Branching ratios (6 13 = 0.22, 8 = rc/2) 




Figure 3: Distribution of the $ ±=l= branching ratios as a function of Aa for #13 = (left panel), 
and #13 = 0.22, 5 = ir/2 (right panel). The asymmetry of the latter plot signals non- vanishing 
CP-violation. 
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3.2.3 Degenerate masses 

When neutrino masses are large compared to the mass differences Am so i and Am atm , all three 
mass states are approximately equal and we can use the model of nearly degenerate neutrino 
masses: mi ~ m 2 ~ m 3 = m. As shown before, the exact value of neutrino mass can not be 
determined by $ ±=t decay statistics, since m becomes independent of branching ratios in the 
degenerate limit. This means, that m is canceled from the expressions for branching ratios and 
the calculation of the Majorana phases is significantly simplified. To obtain general predictions 
we first assume a small value of #13, so that higher order terms of the expansion can be neglected, 
and do not fix other parameters. 

First we can check whether any of the Majorana phases has a non-zero value. When a\ = a 2 = 
0, we would observe nearly equal amount of decays to ee, \i\x and rr channels, while all other 
decay channels would be suppressed: 

BR ee « BR w = BR TT = i (24) 
BR eM = BR eT = BR MT = 0. (25) 

Non-vanishing branching ratios to efi and er channels are clear indicators for non-zero Act. 
When a% = a 2 = a branching ratios to both e/i and er channels are very close to zero BR eM = 
BR eT pa 0. A small non-zero contribution can be added when #13 has a value that is close to its 
upper limit and higher order effects (non-zero sin 2 #13) become influential. 

Very clearly recognizable signature appears when both Majorana phases are maximal («i,2 = tt)- 
If we also assume 823 = 7r/4 (the best fit value), then $ ±=l= has only two possible decay channels 
predicting BR ee = 0.34 and BR^ = 0.66, while all other channels are completely suppressed. 
Small deviations in 823 cause small contributions to the and rr channels while the branching 
ratio to /xr channel is decreased by the same amount. 

The behavior of branching ratios is plotted in Figure H] and Figure [5] which present the de- 
pendence of branching ratios on Aa. The case for a 2 = is shown in the left panel and the 
one for a\ = in the right panel. Figure E] shows the branching ratios for different values of 
Majorana phases when Aa = 0. If we have identified the degeneracy of neutrino masses, we 
can analyze the values of Majorana phases without making any assumption about the values 
of mixing angles. 

• Equal branching ratios to the ee, \x\x and rr channels with all other channels being 
suppressed indicates that a.\ = a 2 = 0. 

• Non-zero branching ratio to the \xr channel means that at least one of the Majorana 
phases has to be non-zero. 

• Non-zero branching ratios to the efi and er channels and the deviation from the result 
BR ee = 0.34 can be generated only by non-zero Aa. Small non-zero contribution to the 
e/i and er channels can be alternatively caused by a large value of #13. 
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Figure 4: Distributions of the $ ±± branching ratios for non-zero Aa. The left figure presents 
the case for fixed a 2 = 0, and the right figure for fixed a\ = 0. 




Figure 5: Distributions of the $ branching ratios for different values of «i and a 2} assuming 
Aa = 0. 



To give exact solutions for the Majorana phases, we fix the values of mixing angles according to 
the tri-bi- maximal model. |Aot| can be found from Eq. (I2U]) which is valid both for the inverted 
mass hierarchy and the degenerate spectrum. Separate values for a,\ and a 2 can be determined 
from the equation system 

_ 2BR Mr - BR ee 1 
Cs = 2BR„„ + BR„ T + BR s „ = 6 (3 - 2C08ai - 4COSQ ' 2) ' 

A • • 4-5C 2 

cosAa = cos ai cos a 2 + sm «! sm a 2 = 77 ttt- (,^o) 

4(1 + C2) 

Again, two possible sets of solutions are found for the Majorana phases and the $ ±=t branching 
ratios do not provide the information to decide which of the solutions is correct. 
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3.3 Effects of non-zero #13 



If our assumption about the exact tri-bi-maximal neutrino mixing should not be valid, the form 
of previously obtained solutions would also be changed. However, the results and methodology 
would generally remain the same. Hopefully new upcoming neutrino oscillation experiments will 
measure the mixing angles with improved precision in the near future jH]. Small changes in 823 
or 612 would not affect the structure of the found solutions and we would only need to substitute 
different values for the mixing angles. Qualitative changes in the analytical expressions appear 
if #13 is taken to be non-vanishing. This influences the structure of the solutions and makes 
the CP-violating angle 6 a physically measurable quantity. We note, however, that, due to 
the smallness of sin #13, the effect of 13 and 5 would enter to the branching ratios as a small 
correction, and extremely precise measurements would be required to detect it. The goal of this 
section is to analyze the effect of non-zero #13 to the previously found solutions for Majorana 
phases and neutrino masses. 

In the following we have still assumed #13 to be a small parameter and considered only the 
leading terms in the expansion with respect to it. We can find the lowest neutrino mass from 
the similar equation as for 13 = (see Eq. (JHJ)), only the measured parameter C[ would differ 
slightly: 

_ 2BR W + 2BR TT + 2Biy - 2BR ee -mj + m\ + 3m§ , 2q 

Cl = 2BR ee + BR e , + BR eT ~ 2m? + m\ + ° (sm ^ (2?) 

Similarly the determination of Majorana phases has exactly the same structure as earlier. Aa 
can be determined uniquely and two possible sets of solutions are found for ai and «2 when we 
attempt to determine the absolute values of Majorana phases. 

In conclusion, assuming small but non-zero #13 does not significantly complicate the determina- 
tion of interesting neutrino parameters at colliders. The solutions would only be slightly more 
complex and involve more decay channels as the relations BR eA1 = BR er and BR^ = BR TT 
no longer hold. There is a theoretical possibility to find solutions also for # 13 and 5, but such 
solutions are very sensitive to experimental errors and, in practice, cannot be used at the LHC. 



3.4 Estimation of the impact of experimental uncertainties 

In this section we consider the effects of experimental uncertainties to the determination of 
$ ±± leptonic branching ratios at colliders and, consequently, to the determination of neutrino 
parameters in collider experiments. The sources of the uncertainties under consideration are. 



• Statistical errors that are relevant for a small number of reconstructed events. In this 
case the number of events observed in particle colliders follow the Poisson statistics with 
theoretically expected average number of events as a mean value. 

• Random measurement errors that dominate in the case of large statistical samples 
and result from the random uncertainties in particle detection and event reconstruction. 
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The measurement errors can vary greatly for different decay channels and their values are 
strongly experiment and detector specific. 

• Systematical measurement errors in particle detection and event reconstruction that 
are negligible for well calibrated detector and correct event reconstruction phenomenology. 



For the numerical simulation of experimental uncertainties we have firstly modified the theo- 
retically expected number of doubly charged Higgs production events N theor with the Poisson 
distribution, then calculated and normalized the corresponding branching ratios and finally 
modified them with Gaussian distortion functions to account for random measurement errors. 
Possible systematical errors have been neglected. Note that for $ ±=l= pair production each 
detected event comprises two doubly charged Higgs decays and N theor = 2N events . 

In reality the measurement errors are different for different decay channels and their values 
depend on the specific detector. As full detector-specific error analysis is out of the scope of this 
paper, we have used uniform uncertainties for all branching ratios for the rough estimation of 
the effect. In particular, we have assumed Gaussian distortion functions with obr = 0.1BR*j eor , 
where BR*^ eor is the theoretically expected branching ratio into the corresponding decay channel 
and i,j = e, //, r. Finally we have run the simulation with randomly distorted branching ratios 
for 50000 times, calculating each time the neutrino parameter of interest. As a result we 
get the distribution function of particular neutrino parameter which measures the stability of 
previously found analytical solutions. 



3.4.1 Mass hierarchy determination 

We remind that the neutrino mass hierarchy is identified by the parameter Ci, defined in 
Eq. dH|), as follows: C\ > 1 corresponds to the normal hierarchy, C\ < 1 to the inverted 
hierarchy and C\ ~ 1 to the nearly degenerate mass spectrum. In general, if the lowest neutrino 
mass is close to zero, the mass hierarchy is very well determined. When the mass increases, the 
distribution of branching ratios is reaching the nearly degenerate limit and the mass hierarchy 
or sign(Am aim ) is increasingly more difficult to measure. 

As an example we have analyzed the behavior of C\ for three different cases: the normal hierar- 
chy for m = 0.02 eV, the inverted hierarchy for m = 0.02 eV and the nearly degenerate limit 
for m = 0.2 eV. The results are presented in Figure El which shows the simulated experimental 
distribution of G\ for two cases with different statistical samples of events. Those imply the 
following la errors. 

• Normal hierarchy (m = 0.02 eV) 

- 1000 $ ±i: decays: d = 6.6 ± 1.1 > 1, 

- 100 $±± decays: C x = 6.6 ± 2.1 > 1. 
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Distribution of parameter Cy 100 decays 
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Figure 6: Simulated distributions of the parameter C\ due to experimental errors for statistical 
samples of 1000 and 100 $ ±± decays in the left and right panel, respectively. We have assumed 
&br — 0.1BR*^ eor for the branching ratio measurement errors. The solid line represents the 



inverted hierarchy with m = 0.02 eV, dot-dashed line the degenerate spectrum with m 
eV and dashed line the normal hierarchy with mo = 0.02 eV. 



0.2 



• Degenerate limit (m = 0.2 eV) 

- 1000 $ ±i: decays: C x = 1.0 ± 0.3, 

- 100 $ ±± decays: C x = 1.0 ± 0.5. 

• Inverted hierarchy (m = 0.02 eV) 

- 1000 $ ±i: decays: C x = 0.06 ± 0.15 < 1, 

- 100 $ ±± decays: C x = 0.06 ± 0.28 < 1. 

The results show that sufficiently good hierarchy detection accuracy is achieved for small lightest 
neutrino masses in both the normal and the inverted hierarchy cases. The accuracy decreases 
when the mass mo increases and C\ — > 1. Such tendency can be understood by comparing the 
normal and the inverted hierarchy plots in Figure [1] where the distribution of branching ratios 
clearly differs for small mass values and becomes very similar to each other when the mass is 
increased. Figure E] shows that the normal hierarchy is very well determined by the parameter 
C\ even for small statistics while at 3a level the inverted hierarchy can be confused with the 
degenerate mass spectrum for small statistics. The main factor that clearly distinguishes the 
normal mass hierarchy with small mo is the negligible value of BR ee that can not occur for the 
inverted or degenerate spectra. 



3.4.2 mo measurement 



After the neutrino mass hierarchy has been determined, we can use either Eq. (1151) (for the 
normal hierarchy) or Eq.( fl6l) (for the inverted hierarchy) to estimate the value of the lowest 
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neutrino mass. In the following we analyze the achievable precision for these parameters. Again, 
for small values of m the measurement precision is sufficiently high and decreases when m 
approaches the degenerate values. 

As already verified by the hierarchy determination accuracy, the normal hierarchy provides a 
distinguishable signature and could thus be easily identified, while the inverted hierarchy can 
be confused with the degenerate spectrum. Such a tendency is also notable in the measurement 
of lowest neutrino mass. For the normal hierarchy not only hierarchy but also the actual value 
of the lowest neutrino mass can be measured with relatively good precision. Turning back to 
our earlier example we assume the true value of the lowest neutrino mass to be mi = 0.02 eV, 
&br — 0.1BR*^ eor to be the branching ratio measurement errors at collider experiments, and 
find that the la experimental errors for the lightest neutrino masses are 5m,i = 2 ■ 10~ 3 eV 
and 8m\ = 5 • 10~ 3 eV for the statistical samples of 1000 and 100 $ ±=t decays, respectively. In 
order to measure from Eq. ( |T6l) for the inverted mass hierarchy with comparable precision, 
a very good statistical basis (more than 5000 events) and the measurement errors smaller than 
&br = O.OIBR-^ 607, are required. The data of such quality would not be obtainable from the 
LHC experiments but only from the future colliders (possibly ILC). 



3.4.3 Measurement of Majorana phases 

Determination of Majorana phases is discussed in detail in Section 13.21 It depends on the 
neutrino mass hierarchy. For the normal mass hierarchy with small mi it is very difficult to 
estimate the Majorana phases with realistic measurement errors. This is due to the fact that 
the only observable «2 does not significantly influence the distribution of branching ratios (see 
Figure [2D- To the contrary, for the degenerate spectrum or inverted hierarchy the Majorana 
phases strongly influence the distribution of branching ratios which, in principle, can be mea- 
sured in realistic experimental conditions. As an example we estimate the measurement error 
for Act for the inverted hierarchy. As we have shown earlier in Eq. (1201) . to high accuracy such 
a calculation does not depend on the value of m 3 . We find that the la errors for Aa are 0.067T 
and 0.037T for 100 and 1000 $ ±± decays, respectively. Similar precision is achieved assuming the 
degenerate spectrum. This result is general and does not depend considerably on the particular 
value of Aa. 

Full detector-specific analysis for the measurement errors of branching ratios requires separate 
analyzes. The error estimations that are found in this section are only approximate, but still 
emphasize the promising nature of our method for determining neutrino parameters in particle 
collider experiments. 



4 Determination of triplet Higgs VEV 

In our scenario the neutrino mass matrix is directly related to the doubly charged Higgs leptonic 
branching fractions according to Eq.(jBJ), and the overall normalization factor is the triplet Higgs 
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VEV v$>. Therefore, one needs additional experimental measurements for determination of v$, 
and thus the entire low energy neutrino mass matrix. Those measurements can come either 
from collider physics, from the low energy neutrino mass measurements or from cosmology 

Let us first assume that v$ is large enough to imply, according to Eq. (JHJ), observable fraction 
of the decays $ ++ — > W + W + , and the collider experiments are sensitive enough to measure 
not just the branching fractions but also the partial widths of the triplet, namely and Tvi/vk- 
The latter may not be possible at LHC but could be possible at ILC experiments [43]. In such 
a case one gets from Eq. (@J and Eq. ffih, 



BRu Tu Tu TuBRww T tot BRww , nc 

= J — — = \ , (2? 



BR^iy Tv^vf kv\ V kBRu V k 

and the determination of v$ from collider experiments is possible. 

If the collider experiments are not able to measure the partial widths of the triplet, one needs 
additional information on the neutrino mass matrix. Assuming that the branching ratio to 
WW channel is measured at any accelerator experiment and |(j7i„) ee | is probed from 0i//3/3 
experiment one gets 

BR ee T ee 

1 \{m v ) ee \ 2 m^±± 



BR^iy Tww 327r kv% 
Now v$ can be directly found as 



i 



[m v ) ee \ z m^±±BR ww ^ ' 



(29) 



V * = { 3^BR^ ) ■ (30) 

Finally, if v $ is too small to imply observable $ — > WW decay rates at colliders, one has to 
rely entirely on leptonic data. If one of the leptonic Yukawa couplings is directly measured in 
the accelerator experiments and |(m„) ee | is probed from the 0u(3(3 experiments, one is able to 
derive the VEV from data. As the simplest example, when T ee is measured, perhaps from the 
resonance at e~e~ collider [44], the VEV can be directly found from 



T " = 32,4 ^" t = V 32-r« • (31) 

As shown, direct measurement of the VEV is possible. However it does require additional 
information which cannot be obtained from the LHC alone. Should the Qv/3/3 yield positive 
results or some of the triplet Yukawa coupling be measured at ILC, we can also give estimates 
on the magnitude of the VEV of the Higgs triplet. 



5 Conclusions 



The main motivation for the present paper is to study how to test the TeV scale triplet Higgs 
neutrino mass mechanism directly at collider experiments. From the collider physics point of 
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view this mechanism has several advantages over the singlet neutrino mass mechanism. As 
the triplet has gauge quantum numbers, its production at colliders is limited only by the mass 
reach not by tiny Yukawa couplings as is the case for singlets. Thus several hundreds of those 
particles can be produced at LHC and ILC experiments. 

The branching ratios of doubly charged Higgs decays to two same charged leptons directly 
probe the corresponding element of the neutrino mass matrix. This allows us to study what 
one can learn about the light neutrino parameters from collider experiments. We have shown 
that the neutrino mass ordering, the lightest neutrino mass and the Majorana phases can be 
measured at colliders by just counting the lepton flavours. We emphasize that those are exactly 
these neutrino parameters which present neutrino oscillation experiments are not sensitive to. 
Therefore collider tests of neutrino mass mechanism may provide a major breakthrough in 
neutrino physics. 

We find that there are some flavour combinations of the doubly charged Higgs decay products 
which definitely point towards certain solutions. For example, should LHC see only doubly 
charged Higgs decays to muons and taus, light neutrinos must have strong normal hierarchy, and 
the lightest neutrino mass can be measured. In particular, the observation or non-observation of 
ee final states is a clear discriminator between the mass hierarchies. Similarly, in the optimistic 
scenarios discussed in Section |HJ one can estimate the magnitude of the Majorana phase(s) of 
light neutrinos. In less clear cases, however, the experimental errors of the collider experiments 
may jeopardize the neutrino parameter measurements and no definite conclusion can be drawn. 

We have also shown that one can actually fully determine the light neutrino mass matrix 
from collider experiments and/or from the measurement of neutrinoless double beta decay 
parameters. This requires determination of the triplet Higgs partial widths to leptons and to 
gauge bosons which could be possible at ILC experiments. If the triplet Higgs turns out to be 
light enough to be produced at colliders, neutrino physics may get an unexpected contribution 
from collider experiments. 



Note added 

When the research presented in this paper was completed, an e-print [15] appeared in arXives 
addressing the same topic. Our numerical results are in agreement with theirs. However, our 
results on neutrino parameters are also obtained in an analytical form which is not the case in 
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